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EHRHART'S POLYNOMIAL FOR EQUILATERAL TRIANGLES IN Z 

EUGEN J. lONASCU 

Abstract. In this paper we calculate the Ehrhart's polynomial associated with a 2-diinensional 
regular polytope (i.e. equilateral triangles) in 1? . The polynomial takes a relatively simple form in 
terms of the coordinates of the vertices of the polytope and it depends heavily on the value d and 



f^ ' its divisors, where d — y ° ^3"^'^ ^nd (a, &, c) (gcd(a, 6, c) = 1) is a vector with integer coordinates 

^Nj ■ normal to the plane containing the triangle. 



■ ■■ 1. INTRODUCTION 

H, 

^ ■ A description of all equilateral triangles with vertices in 1? appeared first in [20j (with the proof 

"^ ' of the full general case in [lOj). An updated version of the same results but with a shorter analysis 

was included in |25j . 

In the 1960's, Eugene Ehrhart (|T^.|15j) proved that given a Z" lattice d-dimensional polytope in 

CN ■ M" (1 < (i < n), denoted here generically by P, there exists a polynomial L('P,t) G Z[t], associated 

> ' 
X/^ ' with V, satisfying 

p 

I> ; (1) L{r, t) = the cardinality of {tV} n Z", te N. 

Lately there has been a lot of literature on various topics generated by the Ehrhart polynomial and 
we cite just a few papers and books: [1], [2], [3], 0], [5], [6], [7], [8], [9], [H], |28], [29], and [3l]. 

r> I Every equilateral triangle in Z^ after a translation by a vector with integer coordinates can be 

^ assumed to have the origin as one of its vertices. Then one can show that the triangle's other 

vertices are contained in a lattice of points of the form (see pOj ) 



(2) Va,b,c '■= {{x,y,z) (^ '^^\ ax + by + cz = 0, a^ + 6^ + c^ = 3(i^, gcd{a,b,c) = 1, a,b,c,d£lj}. 
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Figure 1: The lattice T""^aAc 

This lattice is in general much richer than the sub-lattice, V'^'^a,b,c, of all points which are vertices 
of equilateral triangles with one of the vertices the origin (see Figure [1]) . A more precise description 
of this sub-lattice is given by the following result essentially contained in [25] . 

Theorem 1.1. The sub-lattice V^'^a,b,c is generated by two vectors C ^^'^ V *^ the following 
sense: 7^'^ := AOPQ with P, Q in Va,b,ct is equilateral if and only if for some integers m, n 

(3) O? = m-C -nlf, OQ = n C, +{m-n)lf, with C = (Ci, C2, Cs), 1^ = (?i,?2,?3), "^f = 7, ■> 



(4) 



rac + dbs 



das — bcr i 

C2 = ' ^ 



Xs = r 



ft 



?2 



Zdbr — acs 



3dar + has 



L?3 



where q = a^ + b"^ and r, s can be chosen so that all six numbers in ^ are integers. The sides- 
lengths of AOPQ are equal to d\j2{w? — mn + n^). Moreover, r and s can be selected in such a 
way that the following properties are also true: 

(i) r and s satisfy 2q = s"^ + 3r^ and as a result 2{b'^ + (?) = ^1 + 3Ci and 2(a^ -|- (?) = ^| + 3C| 
(a) r = r'ojx, s = s'ojx where u = gcd{a,b), gcd{r',s') = 1 and x is the product of the prime 
factors of the form 6k — 1 of gcd{d, q) 
(Hi) X divides c 

(iv) \t\ = dy/2, \'t\ = d^/6, and^ -^ = 0. 
(v) s + i\/3r = gcd{A — i^/?,B, 2q), in the ring Z[i\/^], where A = ac and B = bd 
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(vi) This construction is essentially independent of the choice of a and b. In other words, similar 
statements can be formulated for a and c, or b and c instead of a and b. 



Remark: From a computational point of view, it is better to use the smaller two of the numbers 
a, b and c. 

An example here may be illuminating. If d = 15, we observe that we can take a = 1, 6 = 7 and 
c = 25 (3tf2 = a2 + 62 _^ c^). Then, r = -5 and s = 5 give 't = (13, 16, -5) and if = (21, -3,0). 
Also, properties (i), (ii), and (iii) in Theorem 11.11 are satisfied. 

In this paper we are interested in finding the Ehrhard polynomial for this type of triangle in Z^. 
Our formulae seem to be simpler than various other formulae obtained for other polytopes. This 
may be caused by the symmetries and the complications that are already built into our polytopes 
(see the precise formulation below ([S])). 

Prom the general theory (see |3j for a good account) we know that in this case 



L{A,t) = cot^ + yt + C2, te 



where C2 = 1 since we are dealing with a (convex) polytope and cq is the area of the triangle A 
normalized to the area of a fundamental domain of the sub-lattice Va,b,c- We will see that cq is really 
easy to compute but ci, the number of points of the sub-lattice Va,b,c on the sides of the triangle, 
creates some difficulties. Let us point out that this polynomial is not a complete invariant, in the 
sense that we may have the same Ehrhart polynomial for two "different" triangles. By different 
triangles, we understand that one triangle cannot be obtained from the other by the usual trans- 
formations which leave the lattice Z^ invariant. If we take Ai := {(0, 0, 0), (13, —8, 3), (0, —11, 11)} 
and Aa := {(0, 0, 0), (4, 15, -1), (15, 4, -1)}, then 



11 13 

L(Ai, t) = L(A2, t) = —t^ + —t + l,te N. 



The triangles are essentially different since they live in totally different sub-lattices: 7^5,13,13 and 
^1,1,19 respectively. 

It turns out that what we are calculating reduces to counting the number of integer triples 
(x, y, z) satisfying 
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(5) 



' ax + by + cz = 

{3dbr — acs)x — {3dar + bcs)y + qsz > 

[ac{s — 3r) — 3db{r + s)x + [3da(r + s) + bc{s — 3r)]y + (3r — s)qz > 
^ [db{r + s) + ac(s — r)]x — [da{3r — s) + bc(r + s)]y + q{s — r)z + 2qd^t^ > 



,t£n. 



2. A FUNDAMENTAL DOMAIN IN Va,b,c AND THE COEFFICIENT Cq 

In what follows we are going to assume that gcd{a, b,c) = 1 and use the notation introduced in 
Theorem ll.il We notice that from the relations ([3]) and (j4]), we obtain 



(6) ^-^'t-rlf = d{-b,a,0), ^l±il^-Cit = d(0,-c,6), and ^^^^-Cat = d(c, 0, -a). 

It is clear that the vectors it = ^(-6, a,0), it = ■^^^^i-c,0,a) and i^ = ^^rf(5;^(0, -c, 6) 
correspond to points in Va,b,c- We first show that Va,b,c is generated by u, v and w (a point P 
in Pa,6,c is identified by its position vector OP as usual). 

Lemma 2.1. PFii/i the above notation, we have 

'Pa,b,c = {mu + n V -\-pw\m,n,p £ 7j}. 

Proof. If (x, y, z) G Va,b,c-, we see that ax+by+cz = 0. Because uj = gcd{a, b) and gcd{a, 6, c) = 1 
we need to have z = ojz' with z' £ TL. Also, the existence of integers k and I such that ka + ib = to 
is insured by the fact cj = gcd{a, b). This means that we have 

(x, y, z) - z'[gcd{a, c)klf + gcd{b, c)ilj^] = {a,f3',0) e Va,b,c- 

Since aa' + 6/3' = we see that (q',/3',0) = All for some A G Z'^. This shows the inclusion C in 
the equality claimed, and the inclusion D is obvious. ■ 

If we look at the proof of the above lemma we see that it is not necessary to have three vectors 
to generate Va,b,c, but only u and r := gcd{a, c)kv + gcd{b, c)iw , where ka + £b = gcd{a, b), are 
enough. Since there are infinitely many pairs {k, I) satisfying this equality, let us take the solution 
that minimizes k so that k > 0. For computational purposes, we have the following more useful 
result. 
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Lemma 2.2. With the above definition of t we have 

r+s /- r_- 

2d ^ d 



(a) u = ^^ C ~ 5 V 7 where r = ur and s = us (Theorem ] 1.1]) and 



(7) ~^ ^ ~^ T7^' f^^ some a, /3 G Z. 

Proof. The first part follows from Lemma |2. II and (ii) is a consequence of the equalities ^. ■ 
So, u and r form a fundamental domain for Va,b.c- The area of the parallelogram formed by these 
vectors is given by III x r |. 

Lemma 2.3. (i) The area of a fundamental domain for Va.b,c is equal to dv3- 
(ii) The integers a and (3 in ^ satisfy the relation 



(8) \L±liJj^ra\=d 



Proof, (i) We observe that u x v = cdCa c) (^ ^ + ^J + ck) and similarly u 



X vi 



ui gcd(a,c) 

{a i +b j +ck). Hence the area of the parallelogram determined by u and r is equal to 



uj gcd{b,c) 



\lt xl^\ = \-{ak + bl){a't + b~t + c~t)\ = yj a^ + b'^ + c^ = d\^. 

The second part follows from (i) and the equality ([7]), if we take into account that, C, x C, = 
if xlf = 0, C xlf = -if X ( , and | C x ''f I = d'^\/3: 

d^/3 = lit X 7^1 = |S^r X tl = ffi^d^^/3 ^ ®. 

Proposition 2.4. The coefficient cq in the Ehrhart polynomial associated with an equilateral 
triangle 7^'" described by Theorem li.il is given by 

dim? — mn + n^) 
Co = 7, • 



Proof. By the general theory of the Ehrhart polynomial, cq is equal to the area of the trian- 
gle normalized by the area of a fundamental domain of the lattice Va,b,c- Since the area of the 

IdVT 



, ■ 1 ,-r-m,n ■ i , 2d^(m^— mn+n^)\/3 ■ t i o o I i . ■ 2d'^(m^—mn+n'^)^/3 

triangle l^i,\^ is equal to — ^^ -g ^-^^, usmg Lemma 12. 3| we obtam cq — — '^ ^-^^ 



d{m? —mn+n? ) 
2 ■ 
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3. The coefficient ci 




-> 



Figure 2: it ^ ^ ( 



•-:^ 



— > 



^ = fC+§l/ 



We consider the case m = 1 and n = 0. From Lemma l2.lt we see (Figure 2) that the problem 
at hand is to count the points of the lattice, that are determined by u and r , on the sides of the 
triangle T^l^- Clearly we have 3 points, the vertices of the triangle, and in most cases these are 
all of the lattice points on the sides (ci = 3). The points on the sides can be (see Figure 2) on the 
side OP, PQ or OQ. We note that the points on the side OP, are characterized by the existence 
of integers k, i, X and /i such that 



(9) 






Similarly, for points on the side OQ, the equality above changes to 



(10) 



kit + ei^ = {k 






■*^'5'^' 



For the side PQ, we have the characterization 



;ii) 



fc^ + ^^GZ. 



2d 



d 



Lemma 3.1. (i) The coefficient {—k^ +^^) is an integer if £ = fj.^^ + Xr and k = Xfi — fia for 



every A,;U G Z^. 

(a) For the given values of i and k above the other coefficient in ^ becomes 



2d d 

(Hi) //gcd(r, /3) = 1, then there are no points of the lattice in the interior of OP. 
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Proof. Using ([8]), we have 

Similarly, one checks that (ii) is true. 

We observe that ^ implies in particular that gcd(r, /3) divides d. Therefore, if gcd(r, /3) = 1, the 
values in (i) for k and £ give all of the solutions of the Diophantine equation {—k^ + ^f ) = /i € Z. 
(We refer the reader to a basic text on linear Diophantine equations such as |26] or |27j.) Part (ii) 
implies that there are no lattice points on the side OP other than the endpoints. ■ 

Clearly, a similar lemma is true for the sides OQ and PQ. 

So, we have proved the following proposition. 

Proposition 3.2. //gcd(f, /3) = gcd(^,a) = gcd(^,^y^) = 1 then there are no lattice 
points on the sides ofT^'^^ other than the vertices. The Ehrhart polynomial in this case is 

It is natural to ask if all sides of the equilateral triangle 7^ J, ^ may simultaneously contain lattice 
points in their interiors and how many can there be? The next lemma answers these questions. 




Figure 3: All sides may contain lattice points 

Lemma 3.3. (i) If a side contains lattice points in its interior, the intersection of that side with 
T^a,b,c is the set of points which divide that side into d' equal parts, where d' divides d. 

(ii) It is possible that all of the sides of a minimal triangle 7^ ^ ^ have lattice points of Va,b,c 
in their interiors. If di and d2 are the corresponding numbers as in part (i), for two sides, then 
gcd{di,d2) = 1. 

Proof, (i) We refer to Figure 3 in this proof. Without loss of generality, let us assume that 



the closest point to O in Va,b,c on the side OP, is M. Then OM has integer coordinates. Then 
so does kOM for all A; G Z. Therefore, for some k we must have kOM = OP, otherwise we can 
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construct some point in Va,b,ci on the side OP, closer to O than M, using the idea of division with 
non-zero remainder. Therefore, we have shown the first part of (i). For the second part we observe 
that the Diophantine equation ^ = |OMp G Z is possible only if k divides d. This shows (i). 
For (ii), we include here the first example we found of such a triangle, 

A3 = {(0, 0, 0), (220, 539, -539), (747, 12, -267)}, 

corresponding to d = 561 = 3(11)(17), a = 245, h = 613 and c = 713. There are respectively 2, 10 
and 16 lattice points interior to its sides. The associated Ehrhart polynomial is 

,,, , 561*2 + 31t 

L(A3, t) = ^ + 1, t G N. 

For the second part of (ii), let us assume by way of contradiction, that gcd((ii,(i2) = 5 > 1. 
Then there exist integers d'l and d'2 such that di = 6d'i and d2 = 5d'2- From this we get that 

1 d! d! 

1 > ^ = ^ = ^. Without loss of generality we can assume that the side OP is divided into 
di equal parts and the side OQ into d2 equal parts. Then by taking the division points on the 
sides corresponding to d'l and d^ respectively, but different from O, we can construct an equilateral 
triangle in Va,b,c: with side lengths -^ which are strictly smaller than those of 7^ ^ ^. Theorem ll.il 
shows that this is not possible. The contradiction shows that gcd{di,d2) = 1. H 

Let us now look at the case gcd(r, /3) = u > 1. Then Lemma [3T] easily changes into the following. 

Lemma 3.4. (i) //gcd(r, /3) = u, the coefficient {—k^ + ^f ) is an integer if £ = fJ.^-^ + A^ and 
fc = A^ — fia for every A, ^ G Z'^. 

(ii) For the given values of i and k above the other coefficient in (0) becomes 

2d d u 

(Hi) //gcd(r, /3) = u, then there are v — 1 points of the lattice in the interior of the side OP. 

Proof. The calculation we have done in showing Lemma [3. II are valid here with the substitu- 
tion - instead of A. Part (iii) follows from parts (i) and (ii). ■ 
Finally we can put all these things together into the following generalization of Proposition [37 

Theorem 3.5. The Ehrhart polynomial for T^j,^ (given by Theorem \1.1\) is 



'1,0 



^iVbc^t) = TT^ + 1' * e N, where 



dt^ + cit 

+ i 

j2 2 , r 2 , 2 .'- -^ .,^ + 3 , ,/?" — -S 



U^ = a' + ¥ + c\ ci =5cd(r,/3)+gcd(^— ,Q) + gcd(^— ,a + /3), 
and (a, (3) is a particular solution of the Diophantine equation ^-^^/S + ra = d. 
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We observe that this theorem does not depend on the particular solution (a, /3) of the Diophantine 
equation -^/J + ra = d (or equivalently ~^'^^ /S + r{a + /3) = d). 

4. The general case 

We fix two integers m and n with gcd(m, n) = 1 in Theorem ll.il In this case, from ^, we have 
OP = m C, — nif and OQ = n (" + {m — n)~rj . We can solve for C, and 'rf: 

-> {m-n)Op + nO^ , . mO^ - nOp 

C = 2 \ — 2 — °^"'" ^ = ^ \ — 2 ■ 

m^ — mn + n^ m^ — mn + n^ 

Then, the equations ([9]) and ([TO]) change into 

klt + i^ = {k^i + l^)t + (-A;| + £§)^ = 

/ m(r + s) + n(r — s) ma — n(a + /3) \ ^77^ 

(12) Y 2d{m? — mn + n?) d{m? — mn + n^) J 

, n(r + s)-2mr ^ m/3 + na \ 7r± ,77;^ 777^ 
Zdym,'^ — mn + n^ j 0(771^ — mn + n^j / 

Here we have something similar to the case m = 1 and n = 0. 

Lemma 4.1. ('ij T/ie coefficient X in lil!^) is an integer, and equal to t{m — n), if 

k = I3{m?' — mn + n'^)t + [ma — n{a + /3)]i, i = r{m? — mn + n'^)t — (m^^^ + n^^^)t, where t 
and t are arbitrary integers. In this case, the coefficient /x is equal to nt — t. 

(a) The coefficient /i in / f71|) is an integer, and equal to tm, if 

k = —a{m?' — mn + n^)t + (m/3 + na)t, i = ^-^{m? — mn + n^)t + (mr — n^-^)t, where t and t 
are arbitrary integers. In this case, the coefficient A is equal to t — nt. 

(Hi) The value fj, + X is an integer, and equal to tm, if 

k = j3{m? — mn + n^)t + [m(a + /3) — nf3\t, i = r{m? — mn + n^)t — {m^^ + nr)t, where t and t 
are arbitrary integers. In this case, the coefficient fi is equal to nt — t. 

(iv) /f gcd[?7i-^-2^+n^-^,?TT,Q — n(a + /3)] = 1, then there are no points of the lattice in the interior 
ofOQ. 

(v) If gcd{mr — n^-^,m/3 + na) = 1, then there are no points of the lattice in the interior of OP. 

(vi) If gcd{m^Y' +nr,m(a + f3) — n/3) = 1, then there are no points of the lattice in the interior 
ofPQ. 

The proof of this lemma is similar to the proof of Lemma |3. 11 As before, the hypothesis in cases 
(iv), (v), and (vi) in Lemma l4. II can be relaxed. 
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Lemma 4.2. (i) If gcd(r?2^-i^ + ?7,^-^,77ia — n{a + /?)) = v, then the coefficient A is still an 
integer if k = (5{'m? — mn + n'^)t + {ma — n(a + /3))- and i = fire? — mn^r?)t— [m^-^ + n^-^)-, 
where t and t are arbitrary integers. In this case, the coefficient /x is equal to nt— -, and this gives 
V — 1 points of the lattice in the interior of OQ. 

(a) If gcd{mr — n^^^,m(3 + na) = v, then the coefficient fi in \12\l is an integer equal to tm, if 
k = —a{m? — mn + n^)t + (m/3 + na)^, £ = ^^^(m'^ — mn + n^)t + (mr — n^-^^)^, where t and t 
are arbitrary integers. The coefficient A is equal to - — nt. 

(Hi) If gcd(m^Y' + nr,m{a + /3) — n/3) = z^, then the value of fi + X is an integer, and equal to 
tm, if k = j3{m? — mn + n^)t + {m{a + /3) — nfi)-, i = r(rr? — mn + n^)i — {m^^ + nr)-, where t 
and t are arbitrary integers. In this case, the coefficient /x is equal to nt — -. 

(iv) The value v in each of the above cases is always a divisor of d. 

Proof. The calculations we have done in showing Lemma l4.ll are valid here with the substi- 
tution - instead of t. For the part (iv), we note that since gcd(?TT,, n) = 1 we have mx + ny = 1 for 
some integers x and y. Let us consider on (i). One can check the identities 

a I m h n 1 — [ma — n[a + p)) = na, and 



, n\ I ''' + s r - s\ , . ^xx'T-s 

(a + /5) I m— h n— - — 1 + [ma - n[a + p)) —— = md. 

These two equalities can be combined to get 

[[X + y)a + xfi) ( m^ h '^^— ) + ("*" - ^(" + P)) \^^ "^^T' ^ 

which implies the claim that v must divide d. ■ 

The following theorem allows one to compute the Ehrhart polynomial for an equilateral triangle in 
1?. 

Theorem 4.3. The Ehrhart polynomial of an equilateral triangle 7^'" (given by Theorem \l.l\) 
with gcd(m,n) = 1, is given by 

L(X:C.t) = di-^'-mn + n^)t- + c,t ^ ,^ ^ ^ ^^ ^^^^^ 

c\ = gcd(m — - — |-n — - — ,ma—n{a+l3))+gcd{mr—n — - — , m/3+na)+gcd{m — - — \-nr, m{a+(3)—n/3), 
and (a, (3) is a particular solution of the Diophantine equation ^^^13 + ra = d. 
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5. The case of a = 6 and some further questions 

If the equation of the plane dlj) has the property that a = b (2a^+c^ = 3(P), then the parametriza- 
tion of the equilateral triangles given by Theorem 11.11 simplifies to: 

7^'" := AOPQ with P, Q in Va,b,c^ is equilateral if and only if for some integers m, n 

(13) OP = m C - nlf , OQ = nC +im-n)'rf, with C = {(1X2,(3), 'if = ivi,V2,V3), 



(14) 



C2 



d + c 



"■ ^ , and < 



m 



m 



2 

d + c 



We observe that r = s = 1 and so we can choose a = d and /3 = to satisfy ([8|) . Assuming as 
before that gcd(?7i, n) = 1, then the Ehrhart polynomial reduces to the simple formula 



(15) LiTT^'U) 



d{m — mn + n )t + [gcd(?TT., d) + gcd(n, d) + gcd{m — n, d)]t 



+ 1, ten. 



In [23] we have characterized the primitive triples {a,c,d) S N^ satisfying 2a^ + c^ = 3(i^. This 
was done in a manner similar to the way that Pythagorean triples are usually described with a 
one-to-one correspondence to a special set of pairs of natural numbers. 

Theorem 5.1. Suppose that k and i are positive integers with k odd and gcd{k,i) = 1. Then 
a, c and d given by 



(16) d = 2f + k^ with < 



\2f + 2k£ - A;2|, c = \k^ + iki - 2f\, if k ^ i (mod 3) 



a = \2f - 2ki -k^\, c= jfc^ - Aki - 2i\ if k ^ -£ (mod 3) 

constitute a positive primitive solution of 2a^ + c^ = 3d^. 

Conversely, with the exception of the trivial solution a = c = d = 1, every positive primitive 
solution for 2a} + <? = 3d^ appears in the way described above for some I and k. 

In particular, if d > 3 is a prime of the form 8m + 1 or 8m — 5 {m G N), we can find (see 
|13) ) k and i as in Theorem 15.11 So, we can construct examples as in the Introduction in which 
the polynomial is the same for essentially different triangles. As a curiosity, for d = 2011 we get 
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two solutions for {a,c), (913,3235) and (2461,139). The corresponding Ehrhart polynomials are 
identical: LiT^t^t) = ^°^"\+^°^3* + 1, t G N. 

Our ultimate goal is to compute the Ehrhart polynomials for all of the regular lattice polyhedrons 
(tetrahedrons, cubes and octahedrons, see [23]), and naturally the first step is to look at their faces. 

We pose the following question. What is the cardinality of the set 



Eid) := {LiTlll,t)\a^ +b^ + c^= 3d^}7 

This sequence begins in the way recorded in the next table. 

Table 1 



d 


All primitive solutions of a, b, c 


E{d) 


Cl 


1 


{1,1,1]} 


1 


{3} 


3 


{1,1,5]} 


1 


{5} 


5 


{1,5,7]} 


1 


{3} 


7 


{[1,5,11} 


1 


{3} 


9 


{ 1,11, 11], [5, 7, 13} 


2 


{5,11} 


11 


{[1,1,19, 5,13,13, 5,7,17]} 


2 


{3,13} 


13 


{5,11,19, 7,13,17} 


1 


{3} 


15 


{5, 11, 23], [1,7, 25], 5,17,19} 


1 


{5} 


17 


{[7,17,23, 1,5,29, 13,13,23, 11,11,25]} 


2 


{3,19} 


19 


{ 11, 11, 29 , 1, 11, 31 , 5, 23, 23 , 13, 17, 25 } 


2 


{3,21} 


21 


{ 11,19,29, 1,19,31 , 13,23,25} 


1 


{5} 


23 


{ 11, 25, 29 , 1, 25, 31 , 1, 19, 35 , 7, 13, 37 } 


1 


{3} 


25 


{ 1,5,43], [17,25,31], 11,23,35], 5,13,41, 17,19,35]} 


1 


{3} 


27 


{ 1, 31, 35 , 11, 29, 35 , [17, 23, 37 , 7, 17, 43], [13, 13, 43 } 


3 


{5,11,29} 


29 


{ 23, 25, 37 , 1, 11, 49], 1, 29, 41' , 7, 25, 43 , 5, 17, 47 } 


1 


{3} 


31 


{7,25,47, 11,19,49, 5,7,53, 17,35,37, 19,29,41} 


1 


{3} 


33 


{ 5, 29, 49 , 7, 37, 43 , 23, 37, 37 , 19, 35, 41], 25, 31, 41], [23, 23, 47 , [13, 17, 53 } 


3 


{5,15,35} 


35 


{5, 29, 53], [17, 19, 55], 11,23,55, 25,37,41, 25,29,47], 5,13,59} 


1 


{3} 


37 


{7,43,47, 5,19,61], 5,41,49, 1, 25, 59], [11, 31, 55 , 23,37,47} 


1 


{3} 


39 


{ 13, 37, 55 , 1, 29, 61 , 5, 7, 67 , 7, 17, 65 , 11, 31, 59 , 23, 35, 53 } 


1 


{5} 


41 


{ 5, 47, 53 , 1, 1, 71 , 19, 31, 61 , 13, 43, 55 , 25, 47, 47 , [5, 23, 67], 31, 41, 49 , 17, 23, 65 } 


2 


{3,43} 
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